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Abstract. We consider the inverse problem for the wave equation on a 
compact Riemannian manifold or on a bounded domain of IR n , and general- 
ize the concept of domain of influence. We present an efficient minimization 
algorithm to compute the volume of a domain of influence using bound- 
ary measurements and time-reversed boundary measurements. Moreover, we 
show that if the manifold is simple, then the volumes of the domains of influ- 
ence determine the manifold. For a continuous real valued function r on the 
boundary of the manifold, the domain of influence is the set of those points 
on the manifold from which the travel time to some boundary point y is less 
than r(y). 




Figure 1: The grey area depicts the domain of influence and the black curve 
r. The horizontal axis is the boundary of the manifold, and the vertical axis 
represents the direction of the inward pointing normal vector. 
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1 Introduction 



Let M C M. n be an open, bounded and connected set with a smooth boundary, 
and consider the wave equation on M, 

d 2 t u{t, x) - c(x) 2 Au(t, x) = 0, (t, x) G (0, oo) x M, (1) 

u(0, x) = 0, d t u(0, x) = 0, x G M, 

d„it(t, x) = fit, x), (t, x) G (0, oo) x dM, 

where c is a smooth stricty positive function on M, and d v is the normal 
derivative on the boundary dM. 

Denote the solution of Q by u^(t,x) = u(t,x), let T > 0, and define the 
operator 

A2T : / |(o,2T)xaju- (2) 

Operator A 2 t models boundary measurements and is called the Neumann- 
to-Dirichlet operator. Let us assume that c\qm is known but c\m is unknown. 
The inverse problem for the wave equation is to reconstruct the wave speed 
G M, using the operator A 2 t- 
Let T C dM be open, r G C(T), and consider a wave source / in 
L 2 ((0,oo) x dM) satisfying the support condition 

supp(/) C {(t,y) G [0,T]xf; t G [T - r{y), T}}. (3) 

By the finite speed of progation for the wave equation [TBI the solution 
satisfies then the support condition, 

supp(ir (T)) C {x G M; there is y G T such that d(x,y) < r(y)}, (4) 



where d(x,y) is the travel time between points x and see (14) below. Let 
us denote the set in (pi) by M(T, r) and call it the domain of influence. 

The contribution of this paper is twofold. First, we present a method to 
compute the volume of M(T, min(r, T)) using the operator K 2 t- The method 
works even when the wave speed is anisotropic, that is, when the wave speed 
is given by a Riemannian metric tensor g(x) = {gjk{x))^ k , x G M. In the 
case of the isotropic wave equation (JTJ) we have g(x) = (c(x)~ 2 <5jfc)" fc . 

Second, assuming that the Riemannian manifold (M,g) is simple, we 
show that the volumes of M(T, r) for r G C(dM) contain enough information 
to determine the metric tensor g up to a change of coordinates in M. We 
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recall the definition of a simple compact manifold below, see Definition [TJ In 
the case of the isotropic wave equation ([I]) we can determine the wave speed 
c in the Cartesian coordinates of M. 

Our method to compute the volume of the domain of influence is a 
quadratic minimization scheme in L 2 ((0,2T) x dM) for the source / sat- 
isfying the support condition ([3]). After a finite dimensional discretization, 
an approximate minimizer can be computed by solving a positive definite 
system of linear equations. We show that the system can be solved very effi- 
ciently if we use an iterative method, such as the conjugate gradient method, 
and intertwine measurements with computation. In particular, instead of 
solving the equation ([I]) computationally in an iteration step, we measure 
h.2rf for two sources /: one is the approximate minimizer given by the pre- 
vious iteration step and the other is related to the time-reversed version of 



the approximate minimizer, see ( 19 ) below. We believe that our intertwined 
algorithm is more robust against noise than an algorithm where noise is 
propagated by simulation of the wave equation. 

Let us consider next the problem to determine the metric tensor g given 
the volumes of M(dM, r) for all r G C(dM). Our approach exploits the fact 
that C (dM) is a lattice with the natural partial order 

t < a if and only if r(y) < o~(y) for all y G dM. (5) 

Let us define the greatest lower bound of r and a in C (M) as their pointwise 
minimum and denote it by r A a. We recall that a subset of C(M) is a 
meet-semilattice if it is closed under the binary operation A. 
Let us define the boundary distance functions, 

r x : dM^ [0,oo), r x {y) := d(x, y), (6) 

for x G M. We show that the volumes of M(dM, r), r G C(dM), determine 
the meet-semilattice, 

Q(M) = \J{reC(M); r<r x }. (7) 

Moreover, we show that if (M, g) is simple then the boundary distance func- 
tions are the maximal elements of Q(M). The set of boundary distance 
functions determines the Riemannian manifold (M,g) [21, 18J. Thus the 
volumes of M(8M,t), t G C(dM), determine (M,g) if it is simple. 
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Our results give a new uniqueness proof for the inverse problem for the 
wave equation in the case of a simple geometry. Belishev and Kurylev have 
proved the uniqueness even when the geometry is not simple [8]. Their 
proof is based on the boundary control method [2J [7J [T71 [191 121], originally 
developed for the isotropic wave equation [5] . Our uniqueness proof might be 
the first systematic use of lattice structures in the context of inverse boundary 
value problems. 

In previous literature, M(T, r) has been defined in the case of a constant 
function r, see e.g. [18j and the references therein. In this paper we establish 
some properties of M(T, r) when T C DM is open and r G C(T). In partic- 
ular, we show that its boundary is of measure zero. This important detail 
seems to be neglected in previous literature also in the case of a constant 
function r. 

Our method to compute the volume of a domain of influence is related 
to the iterative time-reversal control method by Bingham, Kurylev, Lassas 
and Siltanen Their method produces a certain kind of focused waves, and 
they also prove uniqueness for the inverse problem for the wave equation using 
these waves. Moreover, they give a review of methods that use time-reversed 
measurements [21 HI HH [121 HH US EQ] • A modification of the iterative time- 
reversal control method is presented in [13]. 

2 Main results 

Let (M, g) be a C°°-smooth, compact and connected Riemannian manifold 
of dimension n > 2 with nonempty boundary DM. We consider the wave 
equation 

d$u(t,x) + a(x,D x )u(t,x) = 0, (t, x) G (0, oo) x M, (8) 
«|t=o = 0, d t u\ t=0 = 0, 

b(x, D x )u(t, x) = f(t, x), (t, x) G (0, oo) x dM, 

where a(x, D x ) is a weighted Laplace-Beltrami operator and b(x, D x ) is the 
corresponding normal derivative. In coordinates, {g^ k {x))^ k=1 denotes the 
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inverse of g(x) and \g(x)\ the determinant of g(x). Then 

a(x,D x )u := - ^ ^(ar) -1 |</(a:)|"3^ ^(^Wl V^) J^H 



i,fe=i 

n 

b(x,D x )u : = ^(x)g 3k (x)v k (x 



j,k=l 

where /x is a C°°-smooth strictly positive weight function and 1/ = (z/ 1; . . . , z/„) 
is the exterior co- normal vector of dM normalized with respect to g, that is 
Y^k=i 9^ v 3 v h — 1- The isotropic wave equation (jlj) is a special case of (Tsl) 
with g(x) := (c(x)~ 2 8j k )™ k=1 and ji(x) = c(x) n ~ 2 . 

We denote the indicator function of a set A by 1^, that is, 1a(%) = 1 if 
x G A and 1a (x) = otherwise. Moreover, we denote 



L := {(t, s) G M 2 ; t + s < 2T, s>t> 0}, (9) 
and define the operators 

i r 2T 

Jf(t) := - J l L (t, s)f(s)ds, Rf(t) := f(2T - t), 

K := JA 2T - RA 2T RJ, I fit) : = l (0 , T) (t) / f(s)ds, 

Jo 

where A 2T is the operator defined by being the solution of Q. We 

denote by dS g the Riemannian volume measure of the manifold (dM, g\Q M ). 
Furthermore, we denote by (•, ■) and ||-|| the inner product and the norm of 
L 2 ((0, 2T) x dM; dt <g> dS g ). We study the regularized minimization problem 

argmin ((/, Kf) - 2(1 f, 1) + a \\ff) , (10) 

feS 

where the regularization parameter a is strictly positive and S is a closed 
subspace of L 2 ((0, 2T) x dM). 

Operator a(x, D x ) with the domain H 2 (M) (~)Hq (M) is self-adjoint on the 
space L 2 (M; dV^), where dV^ = nlg^^dx in coordinates. Thus we call dV^ 
the natural measure corresponding to a(x,D x ) and denote it also by m. In 
[9] it is shown that 

(u f (T),u h (T)) LH M; d v,) = (f,Kh). (11) 
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This is a reformulation of the Blagovestchenskii identity [TD] . In Lemma [3] 
we show the following identity 

(v,f(T), l)i2( M . d y M ) = (If, 1). (12) 

This is well known at least in the isotropic case, see e.g. [Bj. The equations 



( [lTp and ggj) imply that 

(f,Kf)-2(If,l) = \\u f (T)\\ 2 L2{M . dVfi) -2(</(T),l) iW v M ) (13) 

= |K( T ) - 1 IL2 (M;d y M) + C, 

where C = — \\l\\if>(M;dv u ) ^ oes no ^ depend on the source /. Thus the mini- 
mization problem ( 10 ) is equivalent with the minimization problem 

argmin (\\u f (T) - l||* 2(M;dVW + a \\ff) . 

For r C M and r : T — > R, we define the domain of influence, 
M(V, r) := {x G M; there is y G F such that d(ar, y) < r(y)}, (14) 



where d is the distance on the Riemannian manifold (M,g). In Section [4] we 
show the following two theorems. 

Theorem 1. Lei a > and let S C L 2 ((0, 2T) x dM) be a closed subspace. 
Denote by P the orthogonal projection 

P : L 2 ((0,2T) x dM) S. 

Then the regularized minimization (iffy has unique minimizer f a G S, and 
f a is the unique f G S solving 

(PKP + a)f = PI + 1, (15) 

where I + is the adjoint of I in L 2 ((0,2T) x dM). Moreover, PKP + a is 
positive definite on S. 

Theorem 2. Let T C dM be open, r G C(T) and define 

S = {f G L 2 ((0,2T) x dM); supp(/) satisfies (g}. 
Let f a , a > 0, be the minimizer in Theorem^ Then in L 2 (M) 

lim u fa (T) = 1 A/(r,rAT) • 
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We denote m T := m{M{dM,r)), for r G C(dM), and := m(M). 
Moreover, we define 

Q(M) := {t G C{dM); m T < m^}, (16) 
R(M) := {r x G C(dM); x G M}, 

where r x is the boundary distance function defined by (|6 |) d be ing the dis- 
tance on the Riemannian manifold (M,g). We denote by Q{M) the closure 
of Q(M) in C(M). In Section [5] we prove the equation §7§ and show the 
following theorem. 

Theorem 3. // (M, g) satisfies the condition 
(G) Xi, x 2 G M and r Xl < r X2 imply Xi = x 2 , 



then R(M) is the set of maximal elements ofQ(M). 

Let T > ma,x{d(x, y)\ x G M, y G dM}. Then the set of volumes, 

V := {m T ; r G C(dM) and < r < T}, (17) 
determines the set Q(M). Note that r x (y) < T, for all x G M and all 



y G dM, and that m-oo = maxV. Moreover, by Theorem [2] and equation (11) 
we can compute the volume m T G V as the limit 

m T = \im(f a ,Kf a ). (18) 

a— w 

The set R(M) determines the manifold (M,g) up to an isometry [2TI ITB"] . 
Hence the volumes (17) contain enough information to determine the man- 



ifold (M,g) in the class of manifolds satisfying (G). In section |5j we show 
that simple manifolds satisfy (G). 

Definition 1. A compact Riemannian manifold (M,g) with boundary is sim- 
ple if it is simply connected, any geodesic has no conjugate points and dM is 
strictly convex with respect to the metric g. 

Let us discuss Theorem [T] from the point of view of practical computa- 
tions. When the subspace S is finite-dimensional, the positive definite system 
of linear equations (15) can be solved using the conjugate gradient method. 
In each iteration step of the conjugate gradient method we must evaluate 
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one matrix- vector product. In our case, the product can be realized by two 
measurements 

A 2 r/, A 2T RJf, (19) 

where / is the approximate solution given by the previous iteration step. 
The remaining computational part of the iteration step consists of a few 
inexpensive vector-vector operations. Thus if we intertwine computation of 
a conjugate gradient steps with measurements (19), the computational cost 
of our method is very low. 



3 The open and the closed domain of influ- 
ence 



Let us recall that the domain of influence M(T, r) is defined in (14) for 
r C M and r : T — y R. We call M(T, r) also the closed domain of influence 
and define the open domain of influence 

M°(r,r) := {x G M; there is y G T s.t. d(x,y) < r(y)}. 

Let us consider the closed domain of influence M(T, r) when T C dM is 
open and r is a constant. Finite speed of propagation for the wave equation 
guarantees that the solution at time T is supported on M(T, r) whenever 
the source / satisfies the support condition Moreover, using Tataru's 
unique continuation result [23, EZ], it is possible to show that the set of 
functions, 

{u f (T); f G L 2 ((0,2T) x dM) and supp(/) satisfies (j3 

is dense in L 2 (M°(T, r)), see e.g. the proof of Theorem 3.16 and the orthog- 
onality argument of Theorem 3.10 in [1 8J . It is easy generalize this for r of 
form 

N 

T(y) = J2 T ArM> y^dM, 

where N 6 N, Tj G R and Tj C dM are open, see [9]. However, the fact that 
M(r, r) \M°(r, r) is of measure zero, seems to go unproven in the literature. 
In this section we show that this is indeed the case even for r G C(T). 

To our knowledge, this can not be proven just by considering the bound- 
aries of the balls B(y, r(y)), y G T. In fact, we give below an example showing 
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that the union of the boundaries dB(y,r(y)), for y G dT, can have positive 
measure. 

Example 1. Let C be the fat Cantor set, M C M 2 be open, g be the Euclidean 
metric, (0,1) x {0} C DM and V = ((0, 1) \ C) x {0}. Then the union 
B := {J ye9T dB(y, 1) has positive measure. 

Proof. The fat Cantor set C is an example of a closed subset of [0, 1] whose 
boundary has positive measure, see e.g. [25]. The map 

$ : (s, a) i-» (s + cos a, s + sin a) 

is a diffeomorphism from R x (0, 7r/2) onto its image in M 2 . The image of 
H := dC x (0, 7r/2) under $ lies in _B. As has positive measure so has 
B. □ 

Lemma 1. Let Y C &e open and Ze£ r G C(T). Then the function 

r r , r (x) := inf (d(x, y) - r(y)) 

is Lipschitz continuous and 

M(r, r) = {xe M; r r , r (x) < 0}, (20) 
M°(r, r) = {x e M; r r>T (x) < 0}. (21) 

In particular, M(T,t) is closed and M°(T,t) is open. 

Proof. Let us define 

r(x) := rr r (x), f(x) := min(d(x,y) — r(y)), 

and show that r = r. Clearly f < r. Let x G M. The minimum in the 
definition of r(x) is attained at a point y G T. We may choose a sequence 
(yj)'jL 1 C T such that ?/j — > yo as j — > oo. Then 

r(x) < d(x, — r(yj) — > f(x) as j — > oo. 

Hence r = r. 

Let us show that f is Lipschitz. Let x G M, and let |/o be as before. Let 
x' G M. Then 

f(x') - f(x) < d(x', y ) - r(y ) - (d(x, y ) - r(y )) < d(x, x). 
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By symmetry with respect to x' and x, r is Lipschitz. 

Let us show (20). Clearly r(x) < for x £ M(T,r). Let x £ M satisfy 
r(x) < 0, and let y$ be as before. Then 

d(x,y ) - r(y ) = r(x) = r(x) < 0, 



and x £ M(r,r). Hence (20) holds. The equation (21) can be proven in a 
similar way. □ 

If T C dM and r is a constant function, then 



M(r,r) = {x £ M; r r , r (x) < 0} = {x £ M; inf d{x,y) < r} 
= {x £ M; < t}. 

Thus for a constant r, our definition of M(T, r) coincides with the definition 
of the domain of influence in [18]. 



Lemma 2. Lei A G M be compact and let t : A — )■ IR 6e continuous. We 
define 

r(x) := inf (d(x,y) — r(y)), x £ M. 

// r zs strictly positive on A or A is a null set, then {x £ M; r(a;) = 0} zs 
a null set. We mean by a null set a set of measure zero with respect to the 
Riemannian volume measure. 

Proof. Denote by V„ the volume measure of M and define 



Let us show that 



Z := {p £ M; r(p) = 0}. 



V g {Z) = V g {Z\A). 



(22) 



If V g (A) = 0, then (|22j is immediate. If r > 0, then 

r(q) < d(q,q) - r(q) = -r(q) <0, q £ A. 



Hence Z n A = and (g holds. 

Let p £ M int . There is a chart ([/, 0) of M int such that <p(p) = and that 
the closure of the open Euclidean unit ball B of lR n is contained in <j)(U). We 
denote U p := <p~ l (B). 
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The sets U p , p G M mt , form an open cover for M int , and as M mt is second 
countable, there is a countable cover U Pj , j = 1, 2, . . . , of M int . Hence 

oo oo 

v g {z) = v g ((z \ A) n (<9M u |J < £ \ A) n 

3=1 3=1 

It is enough to show that 0((Z \ A) D C/ p ) is a null set with respect to the 
Lebesgue measure on B. 

We define for v = (v 1 , . . . , v n ) G M n and x E B, 

n n 
M»(*) '■= Yl vj 9jk( X ) V ^ \ V \ 2 ■= ^(^') 2 , 

j,fc=i j=i 

where (gjk)] t k=i * s the metric (7 in the local coordinates on </>(U). As B is 
compact in 4>(U), there is c p > such that for all t)6R" and x G B 

1 . 

Cp\v\ g(x) < \V\ < —\V\g( x ). 

Cp 

As in the proof of Lemma [T] we see that r is Lipschitz continuous on M. 
Thus by Rademacher's theorem there is a null set N G B such that r is 
differentiable in the local coordinates in B \ N. We denote Z p := <f>((Z\A)n 
U P )\N. 

Let x G Z p and denote p x := As A is compact and q h-> (?) — 

r(g) is continuous, there is g x G A such that 

d{p x , q x ) - r(q x ) = r(p x ) = 0. 

We denote s := d(p x , q x ) = r(q x ). As p x A and q x G A, we have that < s. 

As M is connected and complete as a metric space, Hopf-Rinow theorem 
gives a shortest path 7 : [0, s] — > M parametrized by arclength and joining 
Qx = 7(0) and p x = j(s). For a study of shortest paths on Riemannian man- 
ifolds with boundary see [TJ. As p x G U p , there is a G (0, s) such that 7|[ 0)5 ] 
is a unit speed geodesic of U p C M int . As 7 is parametrized by arclength, 

K7(*)) ^ d (7(*), ?*) - T-(fe) = < - s, * G [a, a], 
and as r(7(s)) = r{p x ) = 0, 

r( 7 (t))-r( 7 (s)) t-a-0 _ 
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The function r o 7 is differentiable at s by the chain rule, and 
^(ro 7 )( g )=lim r(7(t)) - r(7(g)) >l. 

t— >s~ t — S 

As 7 is a unit speed geodesic near s, 

c P = c p \da{s)\ g[x) < \da{s)\ < — \d t j(s)\ g{x) 



Cp Cp 



Hence in the local coordinates in B 



\Dr(x)\ > Dr(x) ■ n^^TI > c P d t{r o j)(s) > c p . 

Let e > 0. There is S(x) > such that 

\r(y) — r(x) — Dr(x) ■ (y — x)\ < c p e\y — x\, y G B(x, S(x)), (23) 

and B(x,5(x)) C B. Here B(x,S) is the open Euclidean ball with center x 
and radius 5. 

The sets B(x,S(x)/5), x G Z p , form an open cover for Z p , and as Z p is 
second countable, there is (xj)^ =1 C Z p such that the sets 

B) :=B( Xj ,5( Xj )/5) 

form an open cover for Z p . By Vitali covering lemma there is an index set 
JcN such that the sets Bj, j G J, are disjoint and the sets 

Bj := B(x j ,5(x j )), j G J, 

form an open cover for Z p . 
We denote 

Dr ( x j) x xi \ 
3 \Dr(xj)y 3 V 3> 



IfyeZ p n Bj, then r(y) = = r(xj) and by (|23j) 

|«i • (y - < \ Dr \ x .^ c P e \y - x i\ ^ e5 r 

We denote by a m the volume of the open Euclidean unit ball in R m and 
by V the Lebesgue measure on B. Let j G J. Using a translation and a 
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rotation we get such coordinates that Xj = and Vj = (1, 0, . . . , 0). In these 
coordinates 



Z p n B 3 C {{y\ y') G R x W 1 " 1 ; \y l \ < e5 v \y'\ < 5,}. 
Hence V{Z P n Bj) < 2e5 j a n ^ 1 5™~ 1 . Particularly, 

V(Z P n Sj) < e— ^(fl,-) = ec n y(^.) ) 
where c n := 2 • 5 n a n _i/o; n . Then 

\ A) n u p )) = v(z p ) = v(z p n |J < £ n 

< ec n ^ y(Sj-) = ec n V(\J B]) < ec n V(B). 

jeJ jeJ 

As e > is arbitrary, V(4>((Z \A)n U p )) = and the claim is proved. □ 

4 Approximately constant wave fields on a 
domain of influence 



Lemma 3. Let f e L 2 ((0,2T) x dM). Then the equation (12) holds. 



Proof. The map h ^ u h (T) is bounded L 2 ((0,2T) x dM) -> L 2 (M), see e.g. 



[22]. Thus is it enough to prove the equation g2j) for / e C c °°((0, 2T) x dM). 
Let us denote 

:= l) L 2(M;dV M )- 

As a(x, -D^) 1 = and b(x, D x )l = 0, we may integrate by parts 

d tV(t) = -(a(x,D x )u f (t), l)i2(M;dV M ) 

= - ((a(x,Z) x )?/(t),i) L 2 (AWfi) - (u f (t),a(x,D x )l) L 2( M . dVfi) ) 
= (b(x, D x )u f (t),l) L 2 (dM . dSg) - (u f (t),b(x,D x )l) L 2 (dM . dSg) 

= (f(t), l)L2{dM;dS g )- 

As div(0) = for j = 0,1, 

v(T) — I / / f(s,x)dS g (x)dsdt 

JO JO JdM 
<-2T r ft 

l(o,r)(*) / f(s, x)dsdS g (x)dt. 



JdM 
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□ 

The proof of Theorem [T] is similar to the proof of the corresponding result 
in [9]. We give the proof for the sake of completeness. 

Proof of Theorem [II We define 

E(f) := (f,Kf)-2(If,l) + a\\ff. 

Then 

E(f) = \\u f (T) - 1||' 2(AW) - ||1|| + « Wft ■ (24) 

Let (fj)f = i C S be such that 

lim E(f,) = mfE(f). 

j^oo feS 

Then 

a || /ill < E(fj) + \\l\\ 2 L 2 (M . dVfi) , 

and (fj)'^ =1 is bounded in S*. As S is a Hilbert space, there is a subsequence 
of (fj)j2. 1 converging weakly in S. Let us denote the limit by e 5 1 and 
the subsequence still by (fj)'jL 1 . 
The map h >-> u h (T) is bounded 

L((0,2T) x 9M) -»■ ff 5/6 " e (M) 

for e > 0, see [22]. Hence ft, h-> u h (T) is a compact operator 

L 2 ((0,2T) x 9Af) L 2 (M), 

and u*i(T) —¥ w°°(T) in L 2 (M) as j — > oo. Moreover, the weak convergence 
implies 

H/ooll < lim infl^ || . 



Hence 



Woo) = lim ||m /j (T) - l||* 2(Jvf . dvy - l|l|li 2(A/Wfi) + a H/ooll 2 



j->oo 



< hm ||^(T) - 1||* 2(JWm) - l|l|| Vw M ) + «l™^f H/,11 2 
= liminf£(/,) = inf£(/), 
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and /oo G S is a minimizer. 

Let / a be a minimizer and h £ S. By orthogonality of the projection P 
and identity (11), it is clear that PKP is self-adjoint and positive semidefi- 
nite. Denote by the Frechet derivative to direction h. As f a = Pf a and 
h = Ph 

= D h E(f a ) = 2(h, PKPf a ) - 2(h, PI+1) + 2a(h, f a ). 
Hence f a satisfies ((TBI). As PKP is positive semidefinite, PKP+a is positive 



definite and solution of (15) is unique. □ 

Lemma 4. Let T C dM be open, r G C(r) and e > 0. T/ien there is a 
simple function 

N 

where N G N, Tj G ffi. and Tj C T are open, swca taat 

r — e < r e almost everywhere on T and 
r t < r on T. 

Proof. As dM is compact, there is a finite set of coordinate charts covering 
dM. Using partition of unity, we see that it is enough to prove the claim in 
the case when T C M n_1 is an open set. But then r is a continuous function 
on a compact set T C lR n_1 , and it is clear that there is a simple function 
with the required properties. □ 

Lemma 5. Let T C dM be open, t G C(T) and let r e , e > 0, satisfy 

t — e < r e almost everywhere on F and 
T f < t on T. 



Then 



limm(M(r,T e )) =m(M(T,r)). 



Proof. Let e > and denote by N C T the set of measure zero where 
r — e > r t as functions on T. Let us show that M°(T, r — e) C M(r, r e ). Let 
x G M°(r, t — e). Then there is ?/o G T such that 

d(x,y ) < r(y ) - e. 
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As r and the function y i— > d(x,y) are continuous and r \ iV is dense in T, 
there is y G T \ N such that 

d(x,y) < r(y) - e < r e (y). 

Hence x G M(T,r t ). A similar argument shows that M(T,r e ) C M°(T, r). 
Clearly M°(r, r - ei) C M°(r, r - e 2 ) for e 1 >e 2 > 0, and 

UM°(r,r-e) = M (r,r). 

e>0 

Hence m(M°(T,T- e)) — >■ m(M°(r, r)) as e — > 0, and 

< m(M°(T, r)) - m(M(r, r e )) < m(M°(T, r)) - m(M°(T, r - e)) 

-> 0, as e -> 0. 

Moreover, by Lemmas [T] and [2] 

m(M(r,r)) = m(M°(r,T)) = limm(M(r, r e )). 



□ 



Proof of Theorem^ We may assume without loss of generality that r < T, 
as we may replace r by r A T in what follows. Let us denote 

S(T,t) := {/ G L 2 ((0,2T) x «9M); supp(/) satisfies Q}. 

By the finite speed of propagation for the wave equation, we have that 
supp(?/(T)) C M(r,r) whenever / G S(T,t). Hence for / G S(T,t), 

||«'(T) - l||* y , = / (^(T) - ifdV, + / W„ (25) 

v ' ^' JM{V,t) JM\M{T,t) 

= \\u f (T) - l M{TjT) \\l 2{M;dV ^+m(M\M(r,T)l 
Let e > 0. By Lemmas [4] and [5] there is a simple function t§ satisfying 
rs<r, m(M(T,T))-m(M(T,T S )) < e. 
By the discussion in the beginning of Section [3j the set 
{u f (T)eL 2 (M(T,Ts)); f G S(T,r 5 )} 
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is dense in L 2 (M(T, t$)). Thus there is f E S(T,Ts) C S(T,t) such that 

\\u (T) - lM(i> 5 )|| L2(JWfi) < e- 



Then 

If 1 1 2 1 1 2 

u (T) - lM(r,r)|| L2(M;dyM) < e + ||lM(r,r s ) - 1 M(r,T)\\ L 2 iM . dV ^ < 2e. 



Moreover, E(f a ) < E(f) and equations (25) and (24) give 
II?/* (T) - l M(r , T )| 



|2 

lL 2 (Af;dV M ) 



= ||u'-(T) - l\\ 2 LHM;dVfi) -m(M\M(T,r)) 

< E(f a ) + ||lfe (M;dVM) -m{M\M{Y,r)) 

< \\uf(T) - 1||* 2(M;(2Vm) - m(M \ M(r, r)) + a ||/|| 2 
= \\u'(T) - lM(r,.)||' 2(AWfi) + a \\f\\ 2 < 2e + a \\f\\ 2 . 

We may choose first small e > and then small a > to get w a (T) arbitrarily 
close to 1a/(i» in L 2 (M). □ 



5 The boundary distance functions as maxi- 
mal elements 

We denote M°(r) := M°(«9M,r), for r E C(dM), and define 

Q:={rGC(M); M\M°(r) ^0}. 

Lemma 6. If t is a maximal element of Q, then r = r x for some x E M . 
Moreover, if the manifold (M,g) satisfies (G), then R(M) is the set of the 
maximal elements of Q. 

Proof. Let x E M and r E C{dM). Then r < r x if and only if x £ M°(r). 
In other words, 

Q = {r E C(dM); there is x E M such that r < r x }. (26) 

Moreover, r x E Q for all x E M. Indeed, r x is continuous and trivially 

r x < r T . 
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Suppose that r is a maximal element of Q. By (26) there is x G M such 
that t < r x , but r x G Q and maximality of r yields r = r x 

Let us now suppose that (M, g) satisfies (G) and show that r x is a maximal 
element of Q. Suppose that rGQ satisfies r x < r. By (26) there is x' G M 
such that t < r' x . Hence 

r x < r < r x ,, 

and (G) yields that x = x' . Thus r = r x and r x is a maximal element of 
Q. □ 

Lemma 7. The set Q is the closure of Q(M) in C(M). 

Proof. Let us first show that Q is closed. Let (t ) )°^ 1 C Q satisfy Tj — > r in 
C(9M) as j ->■ oo. By (g there is (2^)°^ C M such that < r^.. As M is 
compact there is a converging subsequence (xjj^ C Let us denote 

the limit by x, that is, Xj k — > x as k — )• 00. By continuity of the distance 
function, 

r{y) = lim r, fe (y) < lim r Xj (y) = r x (y), y G dM. 

k— »oo fe— >oo K 

Hence t E Q and Q inclosed. 

Clearly Q(M) C Q and it is enough to show that Q(M) is dense in Q. 
Suppose that t e Q. Then there is x G M such that t < r XQ . Let e > 0. As 
M x <9M is compact and the distance function is continuous, there is r > 
such that 

sup \d(x, y) — d(xo,y)\ < e, when d(x,xo) < r and x G M. 
y edM 

Hence r(y) — e < r xo (y) — e < r x (y) for all y G <9M and all x G £>(a;o,r). In 
other words, 

B{x ,r) C (M\M(r-e)), 

and this yields that r — e G Q(M). Functions r — e converge to r in C(dM) 
as e — 7- 0. Thus r is in the closure of A. □ 



Lemmas [6j and [7j together prove Theorem |3j Moreover, (26) yields the 
equation ([7| in the introduction. 

Lemma 8. If (M,g) is simple or the closed half sphere, then (G) holds. 
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Proof. Let X\,x 2 G M satisfy x\ ^ x 2 , and let us show that r Xl ^ r X2 . First, 
if X2 G 9M then r Xl (x2) > = r^^)- Second, if £2 G M mt , then there is 
the unique unit speed geodesic 7 and the unique point y G DM such that 
7(0) = xi, 7(s) = x 2 and 7(5') = y, where < s < s'. As 7 is a shortest 
path from 27 to y (the shortest path if (M,g) is simple) and the shortest 
path from x 2 to y, 

r Xl (y) = s' > s = r X2 (y). 

□ 

The closed half sphere is not simple, since for a point on the boundary 
the corresponding antipodal point is a conjugate point. Hence the manifolds 
satisfying (G) form a strictly larger class than the simple manifolds. 
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